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The concept of a non-equilibrium interfacial tension, defined via the work required to deform 
the system such that the interfacial area is changed while the volume is conserved, is investigated 
theoretically in the context of the relaxation of an initial perturbation of a colloidal fluid towards the 
equilibrium state. The corresponding general formalism is derived for systems with planar symmetry 
and applied to fluid models of colloidal suspensions and polymer solutions. It is shown that the non¬ 
equilibrium interfacial tension is not necessarily positive, that negative non-equilibrium interfacial 
tensions are consistent with strictly positive equilibrium interfacial tensions, and that the sign of 
the interfacial tension can influence the morphology of density perturbations during relaxation. 


I. INTRODUCTION 

The equilibrium interfacial tension between two co¬ 
existing fluid phases is a well established concept of 
enormous relevance to numerous research areas such as 
capillarity, wetting, the morphology of drops and bub¬ 
bles, interface dynamics, and the adsorption of surfactant 
molecules or colloidal particles at interfaces [1, 2]. It is 
well-known that the equilibrium interfacial tension can 
be equivalently viewed either mechanically in terms of 
the difference of the normal and the transversal pressure 
tensor components or thermodynamically via the work 
required to reversibly deform the system in such a way 
that the interfacial area is changed while the system vol¬ 
ume is preserved [3]. Both interpretations require the 
notion of the position of the interface, i.e., of the non- 
uniform part of the system, which can be defined via the 
concept of a Gibbs dividing interface [3]. However, it 
can be shown that the value of the equilibrium interfa¬ 
cial tension between two fluid phases is independent of 
the choice of the Gibbs dividing interface [3]. 

In contrast, the general concept of an interfacial ten¬ 
sion in non-equilibrium situations, e.g., a non-equilibrium 
interface between two coexisting equilibrium bulk phases 
or an interface between two bulk phases being not at co¬ 
existence with each other, has not been uniquely estab¬ 
lished yet. Several attempts have been made to extend 
the notion of an interfacial tension to non-equilibrium 
systems in order to interpret dynamic phenomena at in¬ 
terfaces [4, 5]. The approaches range from time-resolved 
interfacial tension measurements [6] via integrations of 
the Gibbs equation [7-10] to capillary wave analyses 
[11]. The overall assumption underlying all these non¬ 
equilibrium interfacial tensions is that of a sufficiently 
weakly perturbed equilibrium interface. In fact, it has 
been shown recently that the interfacial structure be¬ 
tween two fluid bulk phases not in thermodynamic equi¬ 
librium with each other converges rapidly towards that of 


*Electronic address: bier@is.mpg.de 


the equilibrium interface between the two coexisting bulk 
phases at the same temperature [12]. Hence, different no¬ 
tions of dynamic interfacial tensions can be expected to 
coincide quantitatively provided they do so for equilib¬ 
rium interfaces. However, unlike for equilibrium inter¬ 
faces, the value of the non-equilibrium interfacial tension 
depends on the choice of the Gibbs dividing interface [13]. 

It is well-known that the interfacial tension of equi¬ 
librium fluid-fluid interfaces is strictly positive because 
otherwise the interface would be unstable with respect to 
capillary wave fluctuations. The aim of the present work 
is to demonstrate that this is not necessarily the case for 
non-equilibrium interfacial tensions. To that end, the re¬ 
laxation of an initial non-uniformity of a fluid deep inside 
the one-phase region of the phase diagram is considered. 
Under these conditions the equilibrium state is uniform, 
i.e., no equilibrium interface exists. Therefore, the ini¬ 
tial non-uniformity, i.e., the interface, is a purely non¬ 
equilibrium structure, which ultimately vanishes during 
relaxation towards equilibrium. 

The present analysis is based on the notion of a non¬ 
equilibrium interfacial tension similar to the thermody¬ 
namic definition related to the work of reversible iso- 
choric deformations. It has been shown in Ref. [12] that 
this notion of a non-equilibrium interfacial tension is not 
only consistent with the equilibrium interfacial tension 
but also useful in the context of interfaces between phases 
separated by a fluid-fluid phase transition. In order to 
avoid the difficulty due to the dependence of the non¬ 
equilibrium interfacial tension on the choice of the Gibbs 
dividing interface, setups are considered for which the 
cross-sections of the system perpendicular to one spatial 
direction ( 2 -axis) are all congruent. For such setups the 
number density profiles of the fluid vary only along the 
2 -axis and the interfacial area equals the area of any of 
the congruent cross-sections. Then locating the interface 
position is not necessary. 

In order to apply the definition of the non-equilibrium 
interfacial tension as the work due to a reversible iso- 
choric deformation in practice, the deformation has to be 
performed faster than the relaxation of the fluid under 
consideration. However, this is not possible in practice 
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for simple fluids, since deformations have to be induced 
by deformations of the container walls, from where they 
propagate into the fluid via particle-particle interactions. 
Hence reversible deformations of a simple fluid cannot be 
achieved on time scales faster than the relaxations. 

One way out of this dilemma is to consider fluids of 
colloidal particles which are dispersed in a molecular sol¬ 
vent. In that situation, deformations of the container 
walls induce deformations of the molecular solvent which 
propagate on time scales much shorter than the relax¬ 
ation time of the colloids. Hence the solvent is the tool 
that exerts the force generated by the container walls 
onto the colloidal particles, which are considered as the 
fluid to be studied. 

The general formalism used in this work is derived in 
Sec. H. First some notation to model colloidal suspen¬ 
sions are introduced in Sec. HA. The non-equilibrium 
properties as well as the interfacial structures of col¬ 
loidal fluids can be described simultaneously within dy¬ 
namic density functional theory (DDFT) [14-19], the rel¬ 
evant concepts of which are summarized in Sec. HB. 
The expression of the interfacial tension as a functional 
of the interfacial density profile is derived in detail in 
Sec. HC. The evaluation of the interfacial tension 
for temporally evolving density profiles obtained within 
DDFT leads to the time-dependence of the interfacial 
tension in Sec. HD. The the application of the general 
formalism of Sec. H is illustrated in Sec. HI for three re¬ 
alistic fluid models, all of which exhibit negative values 
of the non-equilibrium interfacial tension. Conclusions 
from the possibility of negative non-equilibrium interfa¬ 
cial tensions and the property of strictly positive equilib¬ 
rium interfacial tensions are drawn in Sec. IV. 


II. GENERAL FORMALISM 
A. Colloidal suspensions 

Consider a three-dimensional monodisperse suspension 
of colloidal particles which interact amongst each other 
via an isotropic pair potential U{r),r = jri — r 2 |. The 
temperature T and the mean number density g of col¬ 
loidal particles are chosen such that in the absence of an 
external field the equilibrium state of the system is that 
of a uniform isotropic fluid well inside a one-phase region 
of the phase diagram. 

The fluid structure of that uniform equilibrium state is 
given by the isotropic direct correlation function c(r) := 
c{r,g),r = jri — r 2 |. It can be shown very generally 
that if t/(r) decays to zero for r —?► oo faster than some 
power law ^ > 5, i.e., if there exists some s > 

5 with lim r^U{r) = 0, the three-dimensional Fourier 

r—>oo 

transform of c(r) is given by 

^q) =Co+%q^ ( 1 ) 

with g = jqj and with TZ{q) decaying to zero for g —>• 0 


faster than g^; the derivative TZ'{q) of TZ{q) with re¬ 
spect to g turns out to decay to zero for g —>• 0 faster 
than g. The situation just described covers practically 
all cases of isotropically interacting colloids, e.g., non- 
retarded dispersion forces in dense suspensions, retarded 
dispersion forces in dilute suspensions, orientationally 
averaged magnetic dipole interaction, as well as pure 
screened Coulomb interaction for suspensions with the 
index of refraction of the solvent being matched to that 
of the colloids. 

At time t < 0 a static external field is applied to the 
colloidal suspension giving rise to a non-uniform equi¬ 
librium number density profile go(r) with mean num¬ 
ber density g = j d^rpo(r). At time t = 0 the 

external held is switched off leaving the colloidal sus¬ 
pension in a non-equilibrium state, which relaxes for 
t ^ oo towards the equilibrium state with the uniform 
number density g. In the following the temporal evo¬ 
lution of the non-equilibrium state of the colloidal sus¬ 
pension is investigated in terms of the time-dependent 
number density prohle p(r,f), which at time f = 0 is 
given by g(r, 0) = £»o(r) and which attains the long time 
limit lim p(r, t) = g. The deviation of the number 

t—^OO 

density p(r, t) from its long time limit g is denoted as 
:= g{r,t) — g, which decays to zero in the long 
time limit: lim 4'{r,t) = 0. 

t—^OO 

B. Dynamic density fnnctional theory 

The relaxation of colloidal suspensions can be de¬ 
scribed as an overdamped conserved dynamics (model 
B), which can be formulated in terms of dynamic den¬ 
sity functional theory [15]. If the system is prepared in 
an arbitrary initial state g{r,t = 0), its state g{r,t > 0) 
evolves with time t such that the Helmholtz free energy 
F[g(t)] reaches a minimum at t —>■ oo. 

The Helmholtz free energy density functional F[g] can 
be written as 

l3F[g] = j d^r g{r) {\n{g{r)A^) - l) + l3F^^[g - g] ( 2 ) 

with the inverse temperature /3 = \/{k'B,T). The excess 
free energy functional /3F®’'[(/)] can be expanded in a func¬ 
tional Taylor series in powers oi (j) = g — g, which results 
in the virial expansion 

j3F^^[(j)\ =,5F®’'[0] — J d^r(j){r) 

+ 0(<(.3). (3) 

The term linear in (j) drops out because J d^r cj){r) = 0, 
and c^‘^\r,g) = c(r, g) = c(r). Since —>• 0 for 
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t —> c», the approximation 

— - yd^r J d^r' c(|r — r'\)(j)(r)(j}{r') (4) 

is applicable, at least at sufficiently late times t. 

Since the colloidal processes considered in the present 
work, which are described in terms of g(r,t), are much 
slower than the molecular degrees of freedom, one can 
assume local thermodynamic equilibrium and define the 
local chemical potential [15] 

Mr,[eW]) = ^bW] 

= ln(p(r,t)A3)+/3;r-(r,[</)(t)]) (5) 

with 

^/r-(r,[0(t)]) = ^[</,(t)] 

Jd^r'c{\r-r'\)(j){r',t). (6) 

The local force —V/r(r, [p(t)]) generates a flux [15] 
j(r, leit)]) = PD (-V^(r, [£>(t)])) g(r,t) 

= -D [<A(t)])) ^.(r, t) 

= -D{Ve{v, t) +QVPg^^{v, [Pit)])) + OiP^) 

=V0(r.i) 

~ -D iVPir, t) + gVPfi^^ir, [Pit)])) (7) 

with the diffusion constant D. From the continuity equa¬ 
tion of the particle number, 

^%^=-V-j(r,[p(t)]), (8) 

one obtains the conserved dynamics (model B [20]) equa¬ 
tion of motion of ^(r, t): 

D {V^Pir, t) + gV^Pn^^ir, [<^(t)])) . (9) 

For the sake of simplicity only planar density profiles 
gir,t) = giz,t), and hence Pir,t) = piz,t), are consid¬ 
ered in the following. Writing 

Piz, *)=■ ^ JPiqz,t) eyiY>iiqzz) (10) 

one obtains from Eq. (6) 

/3/r®^(z, [Pit)]) ~ J dg^c(|g;,|)^(g^,t)exp(%z) 

( 11 ) 


with ciq),q = |q|, denoting the three-dimensional 

Fourier-transform of c(r),r = |r|, (see Eq. (1)). Hence, 
the equation of motion Eq. (9) for piqz,t) reads 

~-Dqlil-gci\qz\))piqz,t), ( 12 ) 

which is readily solved by 

Piqz.t) = ^(9.^,0) exp i-Dqltil - pcdg^l))) 

= ?(<I..0)exp(-^), (13) 


where 


Siq) = 


1 


1 - Qciq) 


(14) 


is the colloidal structure factor of the uniform suspension. 
Equation (13) expresses the relaxation of the Eourier 
mode piqz, t) towards zero on the q^-dependent time scale 


Ti\qz\) 


g(|gz|) 

Dq^z 


(15) 


C. Interfacial tension 

Consider a system with density profiles p(r, t) = 
giz,t), i.e., pir,t) = piz,t), inside a container with all 
cross-sections perpendicular to the z-axis being congru¬ 
ent of area A. The interfacial tension is defined via the 
work which is required to change the interfacial area A 
at constant temperature, volume, and particle number 
[3]. In colloidal suspensions such a change of the interfa¬ 
cial area can be achieved by deforming the incompress¬ 
ible molecular solvent, which drags the colloidal particles 
with it. Due to the separation of molecular and colloidal 
time scales, one is able to perform this deformation, on 
the one hand, sufficiently slowly in order to stay in the 
regime of low Reynolds numbers to avoid dissipation due 
to turbulence, and, on the other hand, sufficiently fast 
such that the colloidal distribution giz,t) is practically 
not evolving during the measurement. 

By slowly deforming the boundaries of the container 
in a time intervall of length tdeform, a flow field v(r, t) 
of the solvent is induced, which fulfills the Stokes equa¬ 
tion within the regime of low Reynolds numbers. Here 
ideform has to be long compared to the relaxation time of 
the molecular solvent (typically nanoseconds) and short 
compared to the relaxation time of the colloids (at least 
milliseconds). Moreover, if the velocity of the container 
walls is negligible compared to the sound velocity of the 
solvent, no pressure gradients occur inside the solvent so 
that the Stokes equation reduces to the Laplace equation 
V^v(r,t) = 0. Complementary to the Eulerian descrip¬ 
tion of the flow in terms of the velocity v(r, t) at position 
r and time t one can use the Lagrangian description in 
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FIG. 1: Sketch of the effect of the displacement field w(r) of 
the molecular solvent induced by a deformation of the con¬ 
tainer walls. Dispersed colloidal particles are dragged along 
with the solvent so that the volume as well as the number of 
contained colloids is conserved for each fluid element (small 
black frame). 


terms of the displacement W (r, t) within the time inter¬ 
val t of a fluid element originally located at position r, 
where both descriptions are related by 

dW 

— (r,t)=v(r + W(r,f),t), W(r, 0) = 0. (16) 

From W(r, t) one obtains the displacement field w(r) := 
W(r, tdeform) experienced by fluid elements during the 
time tdeform- This displacement field w(r) corresponds 
to the deformation map r rw := r-f w(r) of the three- 
dimensional space (see Fig. 1). It is assumed that the 
container is much larger than the colloidal particles so 
that w(r) is slowly varying on colloidal length scales. Due 
to the incompressibility of the solvent, i.e., V- v(r, t) = 0, 

V.w(r)=0(||wf) (17) 

holds for small deformations ||w|| —^ 0. Since the con¬ 
sidered deformations conserve the volume as well as the 
number of colloidal particles of the fluid elements (see 
Fig. 1), the number densities are not changed by the de¬ 
formation map: £iw(rw) = I?(r). 

Upon deformation of (j){r) = g{r) — g with the displace¬ 
ment field w(r), the excess free energy in Eq. (4) leads 
to 





I* W I )0W (I'w )0W w) 

-bw(r) - w(r')|)</>(r)</)(r'). 

(18) 


in Eq. (18) can be expanded for small deformations: 


(;(2)(r — r' + w(r) — w(r')) = c*'^^(r 

-f (wfc(r) - ^^^(r'))^—(r - r') 
OTk 


r') 

0(||wf), (19) 


where an implicit summation over the vector components 
k G {a:, y, z} is performed. Since ^—(r—rO vanishes for 


distances | r — | large compared to colloidal length scales 
whereas w(r) varies smoothly on colloidal length scales, 
the term Wfe(r) —Wk{v') in Eq. (19) can be approximated 
by 


Wk{r) - Wk{r') fv {re - r'^)^ (20) 

with i G {x,y,z}. Therefore Eq. (19) can be rewritten 
as 


(.(2) (j. _ F -t- w(r) — w(r')) = (r — r') 

+ Bke j + ^(ll^ll )’ 

( 21 ) 


where the gradient 


Bke{r) 


dwk 

dri 


(r) 


( 22 ) 


has been introduced. Using the relation 

d 

i{re - r[) exp(iq • (r - r')) = — exp(jq • (r 

oqe 


r')) (23) 


an integration by parts leads to 


dvk 


(r - r'){re - r'g] 




(27r)3 
d3g d 


* 9 fcc(|q|)(rf - r'e) exp(fq • (r - r')) 


^Skec{\q\) -f ^^c'(|q|)^ exp(iq • (r 


(27r)3 


r'))- 

(24) 


Noting i?fcfc(r) = V • w(r) = C>(||w|p) (see Eqs. (17) and 
(22)), one obtains from Eq. (21) 


c(|r — r' -f w(r) — w(r')|) = c(|r — r'|) 



qkqe 


Bke 


/ r -f r' 


Writing c*^^^(r) := c(|r|), the direct correlation function 


c'(|q|)exp(fq- (r - r')) 

+ 0(||wf). (25) 
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Using Eqs. (10) and (25) in Eq. (18) yields 






+ 0(l|wf) 



with i? 22 (p) being the Eourier transform of Bzz{v). 
Whereas - 822 ( 0 , 0 ,^ 2 ) in Eq. (26) contributes only for 
small wave numbers Pz corresponding to those large 
length scales on which the displacement field w(r) varies, 
c'{\<lz\) contributes only for wave numbers Qz correspond¬ 
ing to colloidal length scales. Consequently IP 2 I "C \qz\ 
holds for the dominant contributions to the integrals in 
Eq. (26) so that the approximations 


Therefore, the interfacial tension 7 can be expressed as 
/3 ^[£'w] - PPy] 


B'y = lim 
l|w|H 

47r 


0 Aw — A 
dgz |gz| c'(|g 2 |) ?(gz) 


If ^ ^2 

= -^ h^qzQzc'iqz) (fiqz) 


(31) 


It is worth noting that this expression, according to the 
derivation given above, is independent of the kind of de¬ 
formations of the cross-sections of the container perpen¬ 
dicular to the z-axis, provided the container is larger than 
the size of the colloids such that the displacement field 
induced by a deformation of the container walls varies 
slowly on colloidal length scales. 


D. Time-dependence of the interfacial tension 


^ (=^9z “ y) ~ HFq^z) 


(27) 


apply, thus 


/3F-[0w] = 

+ 0,7*2) J^\qz\c'i\qz\) ^{qz) 

+ 0 (||wf). (28) 


Due to Eqs. (17) and (22), the first integral in Eq. (28), 


r dp 2 

27r 


- 822 ( 0 , 0 ,^ 2 ) 




+ 0(l|wf) 

= -(Aw-A)+ 0 (||wf), 


(29) 


is related to the change Aw — A of the cross-sectional area 
A upon deformation of the container. This renders the 
value of the full free energy Eq. (2) after deformation 

PF[g^] = J d^rw Pw(rw)(ln(pw(rw)A^) - l)+l3F'^^[4>^] 
=Jd^r g{r) (In(p(r)A^) - l) -b PF^^[(j)^] 

= l3F[g] - ——- [dq^ \qz\c'{\qz\) Hqz) 


47r 


+ 0(l|wf). 


Using Eqs. (13)-(15) in Eq. (31) leads to the time- 
dependence of the interfacial tension: 


/37(0 = -y fdqzqzc'iqz) (fiqz^O) 


exp 


2t 

T{qz) 


(32) 


In the following the example of a step-like initial profile 


(j){z,t = 0 ) = 


J -l-Ap /2 , z < 0 
|-Ap/2 ,z> 


0 


(33) 


with Ag > 0 is considered, which leads to the initial 
Eourier transform 


^ Ag 

4>{qzA = 0 ) = i -. 

Qz 


For this case, Eq. (32) reads 

00 

Ag^ 


hit) = — 


27r 


^'( 92 ) 

dg 2 -exp — 

qz 


2t 


iqz) 


(34) 


(35) 


which allows the determination of the time-dependence 
of the interfacial tension hit) by means of the Fourier 
transform c(g) of the direct correlation function. This 
procedure is exemplified in the next section for some fluid 
models. 

The particular shape of the initial profile 4>iz,t = 0), 
i.e., whether it is step-like as in Eq. (33) or smoothly 
varying, can be expected to be irrelevant for the long¬ 
time asymptotic behavior / 37 (t —>■ 00 ). Using Eq. (1) in 
Eq. (35) leads to 


hit) = - 


Ap 2 

27r 


dqz 2 c 2 -f 


n'iqz) 


qz 


exp - 


2t 

hz) 


(36) 


(30) 
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Due to the exponential factor in the integrand, only wave 
numbers q^: with 2t/T{q^) <g; 1 contribute significantly 
to Hence the long-time asymptotics —>• oo) 

is governed by wave numbers q^ corresponding to large 
time scales rijiz). Since inside the one-phase regions of 
the phase diagram 0 < S'(q) < S'(gmax) < oo holds for 
all q € [0,oo), the time scale T{qz) defined in Eq. (15) is 
large only for small wave numbers —>■ 0. For qz ^ 0 
the limiting behavior 


2 c 2 | > 


n'jqz) 

Qz 


and T{qz) 


S{0) 

Dql 


applies, which renders Eq. (36) 


(37) 


/37(t —>■ oo) 


27r 


Aqz 2 c 2 exp ( - 


2Dt 

WY 


-C 2 Ae" 


S{0) _ 

SnDt 


(38) 


Therefore, at long times t —>■ oo, the interfacial tension 
/3j{t oo) is asymptotically proportional to the square 
Ag^ of the density difference Ag and inversely propor¬ 
tional to the square root '/i of the time t. The latter 
property is of course consistent with the ultimate van¬ 
ishing of the interfacial tension due to the disappearance 
of the interface upon relaxation. However, the most in¬ 
teresting feature is that —>■ oo) is proportional to 

the coefficient C 2 of the small-g expansion Eq. (1) of the 
Fourier transform c{q). Depending on the sign of this 
coefficient C 2 , the interfacial tension /37(t oo) can ap¬ 
proach zero with positive or with negative values, i.e., the 
relaxation of the colloidal suspension towards a uniform 
fluid takes place via an interface of low or of high cur¬ 
vature; high curvature interfaces typically exhibit fringes 
(“fingers”). 


III. APPLICATIONS 

In order to illustrate the general formalism derived 
in the previous Sec. H, three examples of model fluids 
are discussed in detail here: polymer solutions, charge- 
stabilized colloids, and colloid-polymer mixtures. 


A. Polymer solutions 


approximation (RPA) c(r) = —j3U(r) one obtains 

'c{qz) = -,5f7o(27r)^/^d^exp 

= ^ + 0{{qzdf)^ , 

(40) 


hence (see Eq. (1)) 

do = -(27r)3/2/3[/od3, (41) 

d2 = V^pUod^. (42) 


Figure 2 displays the temporal evolution of the inter¬ 
facial tension 7 (t) for the case of interaction strength 
/3C/o = 1 and packing fractions rj := ngd^jd S 

{0.05,0.1,0.2,0.5,1}. The thick lines correspond to the 
full expression Eq. (35) based on the direct correlation 
function Eq. (40), whereas the thin lines in Fig. 2(b) are 
the long-time asymptotics Eq. (38). Note that, according 
to Eq. (35), the interfacial tension scales ^ (Arj/r])'^ with 
the initial relative density difference Ag/g = Arj/iq. 

The interfacial tension ')(t) for this system turns out 
to be negative for all times f > 0 (see Fig. 2(a)) because 
c{qz) in Eq. (40) is a monotonically increasing function of 
qz, i.e., c'{qz) > 0, so that, according to Eq. (35), j{t) < 
0. This is, in particular, in accordance with Eq. (38) 
and C 2 > 0 from Eq. (42). Figure 2(b) indicates that 
the asymptotic behavior 'y{t) ^ jjj gq_ ^ 33 ^ 

expected to apply at times beyond the Brownian time 
scale t > d^/D; this statement applies to a wide range of 
the parameters fdUo and 77 , the results of which are not 
shown here. 

According to Refs. [23, 24], the interaction poten¬ 
tial Eq. (39) with interaction strength fdUo = 1.58 and 
d = Rg y^2/3 = 60 nm is valid to describe the interaction 
of pairs of poly(a-methyl styrene) chains with molecular 
mass Mu, = 2.96 x 10® g/mol in toluene at 25 °C in the di¬ 
lute regime rj < 0.22. At the Brownian time t = d^/D = 
46 ms, where D = 0.079/rm^/s is the experimental dif¬ 
fusion constant, one finds (i^{t)d^ = —Q.li{Ar]/'q)'^ for 
77 = 0.22, i.e., for an initial perturbation of IA 77 I = O.I 77 
the free energy k-^T is gained by interface fluctuations 
which increase the interfacial area by 2.77 . 


As the most simple example consider a polymer solu¬ 
tion described within the Gaussian core model [21], i.e., 
two polymer chains interact via the Flory-Krigbaum po¬ 
tential [22] 

/3H(r) = /3Hoexp(^-^^ (39) 

with [/q > 0 measuring the interaction of two polymer 
chains located at the same position and d correspond¬ 
ing to the radius of gyration. Within the random-phase 


B. Charge-stabilized colloids 

As a second example consider a colloidal suspension 
of charge-stabilized colloids of diameter d in a solvent 
with inverse Debye length k. Here the interaction of two 
colloidal particles is modeled by the potential 

{ 00 , r < d 

cirr exp(-K(r - d)) (43) 

- TJd - 
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reference system is given by 
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FIG. 2: Non-equilibrium interfacial tension 7 (t) during the 
relaxation of polymer solutions described by a Gaussian core 
model Eq. (39) with polymer number densities g correspond¬ 
ing to packing fractions t] = TVQcfi/6 € {0.05,0.1,0.2,0.5,1}. 
Thick lines represent the full expression Eq. (35), whereas 
thin lines in panel (b) correspond to the long-time asymp¬ 
totics Eq. (38). Panel (a) shows that, for this system, the 
interfacial tension is negative, 7 (t) < 0, for all times t > 0. 
From panel (b) one can infer that the asymptotic expression 
Eq. (38) applies for times t > /D. 


c(r) 


CHS(?') 

exp(-K(r- d)) 

-pUd - -n - 

r/d 


,r < d 
,r >d, 


(44) 


where the Percus-Yevick direct correlation function of 
hard spheres with diameter d and packing fraction rj is 
[25] 


CHS(?') 


with 


Chs(^) + Chs(^)^ + Chs(^) (^) ^ 

0 ,r > d 

(45) 


Chs(^) 

Chs(^) 

Chs(^) 


(l + 2^)2 

{1-^r ’ 

6?7(1 -I- ri/2)'^ 

??(1 + 2 ??)^ 
2(l-r,)4 ■ 


(46) 


Figure 3 displays the temporal evolution of the inter¬ 
facial tension ^{t) for the case of interaction strength 
PUd = 1, inverse Debye length nd = 1, and packing 
fractions r] G (0.1,0.2,0.3,0.4,0.5}. The thick lines cor¬ 
respond to the full expression Eq. (35) based on the di¬ 
rect correlation function Eqs. (44)-(46), whereas the thin 
lines in Fig. 3(b) are the long-time asymptotics Eq. (38). 

The interfacial tension ^{t) for charge-stabilized col¬ 
loids, as for polymer solutions in the previous Sub¬ 
sec. Ill A, is negative, 7 (t) < 0, for all times t > 0 
(see Fig. 3(a)), which is also in agreement with the posi¬ 
tive values C 2 > 0 entering in the long-time asymptotics 
Eq. (38). Figure 3(b) indicates that the asymptotic be¬ 
havior jit) ^ in Eq. (38) sets in, as for polymer 

solutions in the previous Subsec. Ill A, at times beyond 
the Brownian time scale t > d^/D. 

It has to be noted that the interaction strength pUd 
of some real charge-stabilized colloids can exceed unity 
by orders of magnitudes: In Ref. [26], e.g., an aqueous 
dispersion = I6I nm) of polystyrene sulfate spheres 
of diameter d = 652 nm and packing fraction </> Ri 0.04 is 
reported with pUd = 463. For such a large value of pUd 
the RPA Eq. (44) is not justified so that one has to use 
more reliable schemes to evaluate c(r). 


C. Colloid-polymer mixtures 


where the strength pUd > 0 of the double layer repulsion 
increases with the charge of the colloids. The direct cor¬ 
relation function of the colloidal suspension within RPA 
of the double layer repulsion in excess to the hard-sphere 


As a final example consider a colloid-polymer mix¬ 
ture composed of colloidal hard spheres of diameter d 
suspended in a solution of polymer coils of diameter 
i. Within the Asakura-Oosawa-Vrij theory [27-29], the 
polymer coils with packing fraction r]p give rise to an 
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FIG. 3: Non-equilibrium interfacial tension 7 ( 1 ) during the 
relaxation of suspensions of charge-stabilized colloids inter¬ 
acting via the potential Eq. (43) with strength fiUd = 1 of 
the double layer repulsion, inverse Debye length nd = 1, and 
packing fractions rj £ {0.1, 0.2, 0.3, 0.4, 0.5}. Thick lines rep¬ 
resent the full expression Eq. (35), whereas thin lines in panel 
(b) correspond to the long-time asymptotics Eq. (38). Panel 
(a) shows that, for this system, the interfacial tension is neg¬ 
ative, 7 (t) < 0, for all times t > 0. From panel (b) one can 
infer that the asymptotic expression Eq. (38) applies for times 
t > d?ID. 

effective depletion interaction 

(47) 

between colloidal particles for distances in the range r £ 
(d, d +1). The direct correlation function of the colloidal 


suspension within RPA of the depletion interaction in 
excess to the hard-sphere reference system is given by 

{ CHs(r) ,r<d 

-/3P(r) ,rG(d,d-b£) (48) 

0 ,r >d + 

Figure 4 displays the temporal evolution of the interfa¬ 
cial tension 7 (t) for the colloidal packing fraction rj = 0 . 2 , 
the size ratio ^|d = 2 , and the polymer packing fractions 
ryp G {0,0.01,0.02,0.05,0.1,0.2}. The thick lines corre¬ 
spond to the full expression Eq. (35) based on the direct 
correlation function Eq. (48), whereas the thin lines in 
Fig. 4(b) are the long-time asymptotics Eq. (38). 

In accordance with Eq. (38), the asymptotic decay 
7 (t) —>■ 0 is from below for weak depletion attractions, 
r/p G {0,0.01,0.02}, with positive values C 2 > 0, whereas 
it is from above for strong depletion attractions, rj^ G 
{0.05, 0.1, 0.2}, with negative values C 2 < 0. This demon¬ 
strates that it can be possible to alter the sign of the 
asymptotic decay by simply changing the polymer con¬ 
centration. Moreover, the cases r]p G {0.05, 0.1} in Eig. 4 
show that the sign of the interfacial tension 7 ( 1 ) < 0 at 
early times t —>■ 0 has no influence on the sign j{t) > 0 
at late times t —>■ 00 . Figure 4(b) indicates that also in 
these systems the asymptotic behavior ^{t) ^ in 

Eq. (38) can be expected to apply at times beyond the 
Brownian time scale t > d^/D. 

IV. CONCLUSIONS 

It has been shown in Sec. II that, under the mild condi¬ 
tions of an interaction potential U (r) decaying for r —> 00 
faster than ^ r® with s > 5 (see Subsec. II A) and of 
an underlying model B dynamics (see Subsec. IIB), an 
initially perturbed colloidal fluid relaxes towards equilib¬ 
rium accompanied by a time-dependent non-equilibrium 
interfacial tension 7 ( 1 ) which is not necessarily positive. 

Here the non-equilibrium interfacial tension ^(t) is de¬ 
fined via the work required to deform the system such 
that the interfacial area changes while the fluid volume 
is preserved (see Subsec. IIC). An experimental deter¬ 
mination of the non-equilibrium interfacial tension 7 (f) 
requires, on the one hand, a sufficiently slow dynamics 
of the system such that the structure does not change 
appreciably during the measurement and, on the other 
hand, a sufficiently fast dynamics to preserve local equi¬ 
librium during deformation. Colloidal suspensions and 
polymer solutions with their wide separation of solvent 
and solute time scales (nanoseconds vs. milliseconds) are 
systems for which the negative non-equilibrium interfa¬ 
cial tension could be realized. There the fluid, whose 
non-equilibrium interfacial tension is to be measured, is 
formed by the colloidal particles, whereas the solvent acts 
as a medium whose Stokes flow due to agitation of the 
container walls leads to the deformation of the colloidal 
fluid. 
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FIG. 4: Non-equilibrium interfacial tension 'y{t) during the 
relaxation of colloid-polymer mixtures described within the 
Asakura-Oosawa-Vrij model Eq. (47) with colloidal packing 
fraction rj = 0 . 2 , polymer-to-colloid size ratio £/d = 2 and 
polymer packing fractions rjp € {0,0.01,0.02,0.05,0.1,0.2}. 
Thick lines represent the full expression Eq. (35), whereas the 
thin lines in panel (b) correspond to the long-time asymp¬ 
totics Eq. (38). Panel (a) demonstrates that the interfa¬ 
cial tension 'y{t) as a function of time t can be positive 
( 7 ( 1 ) > 0 , see case r]p = 0 . 2 ) or negative ( 7 ( 1 ) < 0 , see 
cases rjp £ { 0 , 0 . 01 , 0 . 02 }) or it can change its sign (see cases 
rjp £ {0.05,0.1}). From panel (b) one can infer that the 
asymptotic expression Eq. (38) applies for times t > d^/D. 


The non-equilibrium interfacial tension 7 (t) is, under 
the above-mentioned conditions, determined by the equi¬ 
librium structure, the diffusion constant, and the initial 
density difference Ag from equilibrium (see Subsec. IID). 
In particular, Eq. (35) leads to 7 ( 1 ) ~ Ag"^, which is 
reminiscent of the expression for the interfacial tension 


of equilibrium interfaces within the square-gradient ap¬ 
proximation [3]. However, the square-gradient approxi¬ 
mation is applicable only, if the square-gradient contri¬ 
bution, and hence the interfacial tension, is positive, be¬ 
cause otherwise all uniform bulk states would be unstable 
with respect to density fluctuations. A positive square- 
gradient contribution corresponds to a negative coeffi¬ 
cient C 2 < 0 of the quadratic term ^ in Eq. (1), but 
this term is not required to be negative beyond square- 
gradient approaches, where the higher-order terms 7^(g) 
are not neglected. The interesting and rather general 
finding Eq. (38) shows that the long-time limit 00 ) 

is proportional to this coefficient C 2 of the quadratic con¬ 
tribution ^ in Eq. (1). The reason is that the re¬ 
laxation time T{q) in Eq. (15) decreases with increasing 
wave number q so that the contributions of large wave 
numbers q decay quickly (see Eq. (35)). 

Several examples of realistic systems have been pro¬ 
posed in Sec. Ill for which negative non-equilibrium inter¬ 
facial tensions ■j(t) during relaxation can be expected to 
occur: Polymer solutions (Subsec. Ill A, Eig. 2), charge- 
stabilized colloids (Subsec. IIIB, Fig. 3), and colloid- 
polymer mixtures (Subsec. IIIC, Fig. 4). The latter 
system even offers the possibility to switch between an 
asymptotically positive {jit —>■ 00) > 0 ) and an asymp¬ 
totically negative (7(1 —?► 00) < 0 ) non-equilibrium inter¬ 
facial tension. As a rule of thumb, C 2 < 0, and hence 
7(t —>■ 00) > 0 , corresponds to an interaction poten¬ 
tial U{r) with the attractive contribution being suffi¬ 
ciently strong as compared to the repulsive contribution, 
whereas C2 > 0 , and hence j{t —>• 00) < 0 , is the re¬ 
sult of an interaction potential U(r) with the repulsive 
contribution being sufficiently strong as compared to the 
attractive contribution. 

Positive and negative values of the non-equilibrium in¬ 
terfacial tension 'y(t) lead to different morphologies of 
density perturbations: The relaxation towards a uniform 
equilibrium state takes place by minimizing the interfa¬ 
cial area for 7 ( 1 ) > 0 and by maximizing the interfacial 
area for j(t) < 0. Hence, localized density perturbations 
tend to spatially shrink with interfaces being smooth for 
7 (t) > 0 while they tend to spread out with increasingly 
rough interfaces for 7 ( 1 ) < 0 . 

It is important to realize that the relaxation of a den¬ 
sity perturbation is not driven by the non-equilibrium 
interfacial tension, which can be positive or negative 
(Fig. 4(a)), but by the non-uniformity of the local chem¬ 
ical potential (see Eq. (5)); the interfacial tension is 
merely related to the interfacial structure formed dur¬ 
ing the relaxation process. It is straightforward to show 
with Eqs. (5), (7), and ( 8 ) that the Helmholtz free energy 
F[g(t)] is a monotonically decreasing function of time t, 


mm] 

dt 


1 

Jd 


d^r 


JM! 

g{r,t) 


< 0 , 


(49) 


and it is this decrease which leads to the irreversibility 
of the relaxation towards equilibrium. In contrast, the 
non-equilibrium interfacial tension in Eq. (31) quantifies 
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the linear response of the system to reversible deforma¬ 
tions, and it has been shown in Sec. IIC that this linear 
response is independent of the type of deformation. Be¬ 
yond linear response one may find a dependence of the 
system’s response on the type of deformation, but such a 
behavior has to be described by a different quantity than 
the interfacial tension. 

Based on these comments, one can reconcile the non¬ 
equilibrium interfaces during relaxation studied in the 
present work with equilibrium interfaces between two 
coexisting bulk phases: First, equilibrium interfaces, 
whose interfacial tensions do not vanish, exhibit a time- 
independent structure because the local chemical poten¬ 
tial in equilibrium is uniform. Second, bulk coexistence is 
possible only for a sufficiently strong attractive contribu¬ 
tion to the interaction potential U{r) as compared to the 
repulsive contribution, which typically leads to C 2 < 0, 
i.e., 7 (t —?► oo) > 0. Conversely, a negative equilibrium 
interfacial tension requires an interaction potential U (r) 
with the repulsive contribution dominating the attractive 
one, but under these condition no bulk coexistence, and 
thus no equilibrium interface, occurs. 

In conclusion, a rather general expression for the in¬ 


terfacial tension of a relaxing planar non-equilibrium in¬ 
terface is derived and applied to realistic systems of col¬ 
loidal suspensions and polymer solutions. It is shown 
that the non-equilibrium interfacial tension is not nec¬ 
essarily positive, that negative non-equilibrium interfa¬ 
cial tensions are consistent with strictly positive equi¬ 
librium interfacial tensions, and that the sign of the in¬ 
terfacial tension can influence the morphology of den¬ 
sity perturbations during relaxation. The present study 
highlights that the useful concept of a non-equilibrium in¬ 
terfacial tension shares some but not all properties with 
the equilibrium interfacial tension. Until now concepts of 
non-equilibrium interfacial tensions have been introduced 
only for systems close to equilibrium, whose known relax¬ 
ation dynamics towards equilibrium allows for generaliza¬ 
tions of the equilibrium interfacial tension. It is a future 
task of enormous relevance to understand the properties 
of non-equilibrium interfaces also far away from equilib¬ 
rium. Whether the interfacial tension is a useful concept 
also far away from equilibrium or whether its applica¬ 
bility is restricted to the vicinity of equilibrium is an 
interesting open question. 
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